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Abstract: We analyse the process of reheating the Universe in the electroweak theory 
where the Higgs field plays a role of the inflaton. We estimate the maximal temperature of 
the Universe and fix the initial conditions for radiation-dominated phase of the Universe 
expansion in the framework of the Standard Model (SM) and of the i/MSM — the minimal 
extension of the SM by three right-handed singlet fermions. We show that the inflationary 
epoch is followed by a matter dominated stage related to the Higgs field oscillations. We 
investigate the energy transfer from Higgs-inflaton to the SM particles and show that the 
radiation dominated phase of the Universe expansion starts at temperature ~ (3 — 15) x 
10^^ GeV, where the upper bound depends on the Higgs boson mass. We estimate the 
production rate of singlet fermions at preheating and find that their concentrations at 
are negligibly small. This suggests that the sterile neutrino Dark Matter (DM) production 
and baryogenesis in the z^MSM with Higgs-driven inflation are low energy phenomena, 
having nothing to do with inflation. We study then a modification of the i/MSM, adding 
to its Lagrangian higher dimensional operators suppressed by the Planck scale. The role of 
these operators in Higgs-driven inflation is clarified. We find that these operators do not 
contribute to the production of Warm Dark Matter (WDM) and to baryogenesis. We also 
demonstrate that the sterile neutrino with mass exceeding 100 keV (a Cold Dark Matter 
(CDM) candidate) can be created during the reheating stage of the Universe in necessary 
amounts. We argue that the mass of DM sterile neutrino should not exceed few MeV in 
order not to overdose the Universe. 



Keywords: inflation, physics of the early universe, dark matter, cosmological neutrinos| . 



Contents 

|l]. Introduction 

|2|. Higgs-driven inflation 

|3|. Reheating in Higgs-driven inflation 
|3.1| Qualitative picture 
|3.2| W boson production 



3.3 Transfer into relativistic particles 



3^ Higgs production 
5| Fermion production from Higgs decay 



^ Initial conditions for the hot Big Bang 

|5|. Higher dimensional operators 

[l| Contributions of the higher dimensional operators to inflation 
Higgs operators 



5.1.1 



5.1.2 Yukawa terms 



5.2 Sterile neutrino production 



5.2.1 Thermal production 



5.2.2 Production during preheating 



|6|. Higher dimensional operators and baryon asymmetry 

^ Conclusions 

|A| . W boson production 
[A.l Semiclassic approach 



A. 2 Non-resonance production 



A. 3 Stochastic resonance 



W boson decays and scatterings: energy transfer into relativistic parti- 
cles 



12 
13 

n 

15 
16 

m 
m 

m 



24 
25 
26 

Be 



B.l| W boson decays 

p3.2| W bosons annihilation 

1^. Non-resonant Higgs production on the nonlinearities of the potential for 
small X 

p| . Tunnelling through a — |x| barrier 
^ Useful sums 

1^. Nonresonant particle production with h and \h\ sources 



30 

m 



32 
33 
34 
35 



- 1 - 



1. Introduction 

The statement that the Universe was dense and hot in the past is an estabUshed exper- 
imental fact. It follows from existence of the Cosmic Microwave Background radiation 
(CMB), which has a perfect Planck spectrum, and from accordance of predictions of the 
Big Bang Nucleosynthesis (BBN) with observations. The latter tells that the Universe 
had the temperature of at least few MeV. Whether the Universe was even hotter is an 
open question, which hardly can be answered by experimental means and thus is biased 
by theoretical prejudice. 

According to the current theory of inflation (for a recent review see ) the early evo- 
lution of the Universe can be roughly divided into three parts. During the first, inflationary 
stage, the Universe expands exponentially and becomes nearly flat. At this stage relic grav- 
ity waves and matter perturbations, leading to structure formation, are generated. During 
the second, reheating stage, the energy stored in the inflaton field is transferred to the 
fields of the Standard Model and other (hypothetical) particles (if they exist). The third 
stage is the radiation dominated Universe in nearly thermal equilibrium for most of the 
SM particles. The starting moment of this stage tr corresponds to a maximal temperature 
of the Universe Tmax) and this is the onset of the standard Hot Big Bang. 

The system in thermal equilibrium is completely characterised by temperature T and 
chemical potentials for exactly conserved quantum numbers Qi; the corresponding op- 
erators Qi obey [Qj,?Y] = 0, where Ti. is the Hamiltonian of the system. For the expanding 
Universe the precise thermal equilibrium never exists. To describe the state of the Universe 
at T the set of operators Qi should be supplemented by approximately conserved 

operators Qa, whose rate of change is much smaller than the rate of the Universe ex- 
pansion. Thus, to follow the Universe evolution at later times, t > tr, one can use the 
ordinary kinetic approach based on Boltzmann equations (or equations for density matrix, 
if coherent quantum effects are essential) with initial density matrix 



The magnitude of the maximal temperature Tmax together with the set of values of the 
chemical potentials ;Uj, can be called the initial conditions for the Hot Big Bang. If they 
are known, the further evolution can be completely specified by the standard methods of 
kinetic theory. 

Clearly, to find the initial conditions for the Big Bang one has to know what are the 
relevant particle degrees of freedom at T < Tmax (in particular, if any new particles beyond 
those already present in the SM exist), or, in other words, what is the Hamiltonian 7i. The 
knowledge of the Hamiltonian would allow to determine the set of conserved Qi or nearly 
conserved Qa operators and identify the relevant chemical potentials. Now, to determine 
Tjnax and ^ii^A the interaction of the inflaton with the fields in 7i must be known, and the 
physics of reheating must be elucidated. 

Basically, to find the initial conditions for the Big Bang one should have at hand the 
theory which is valid up to the scale of inflation. There are quite a number of proposals 
for these types of theories, based on different ideas about physics beyond the SM. These 
ideas include low energy supersymmetry and Grand Unification, small, large or infinite 
extra dimensions (see e.g. and Q for reviews) and many others. Clearly, any model of 




(1.1) 
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physics beyond the SM, must be able to explain the observed phenomena that cannot be 
addressed by the SM physics. They include neutrino masses and oscillations, the existence 
of dark matter in the Universe, baryon asymmetry, inflation, and accelerated expansion of 
the Universe at present.^ 

The most economical particle physics model which is capable of solving in a unified 
way all these problems of the SM is the i/MSM (Neutrino Minimal Standard Model) of 
1^, ^. This theory is nothing but the SM augmented by three relatively light (lighter 
than Z boson) right-handed singlet fermions. If the dilaton field is added to the i^MSM, 
the theory can be made scale-invariant at the quantum level by a specific renormalization 
procedure ||6|, |^. The spontaneous breaking of the scale invariance leads then to generation 
of all mass parameters, including the Newton's gravity constant. Higgs mass is stable 
against quantum corrections, cosmological constant is equal to zero, while dark energy, 
leading to the late acceleration of the Universe, appears if general relativity is replaced by 
the unimodular gravity [^. Different phenomenological and cosmological aspects of this 
theory, together with the study of how to search for new particles, can be found in refs. 
@, |, |, 0, |n|, |1|, [l^, |1|]. In it was argued that this model may 



be valid all the way up to the Planck scale (for a similar argument in a related theory, 
see ref. [pOt]). Many of the parameters of this model are already fixed or constrained by 
existing cosmological observations and particle physics experiments. 



In 1 21 1 it was found that the Higgs boson of the SM can play a role of the inflaton, 
if its non-minimal coupling to the gravity Ricci scalar is large enough. Exactly the same 
mechanism works in z^MSM and its scale-invariant version with dilaton The evolution 
of the dilaton in the latter model with phenomenologically interesting choice of parameters 
happens to be irrelevant for inflation. 

Reference [^] provides a rough upper limit on the maximal temperature of the Uni- 
verse. The aim of this work is to demonstrate that the problem of the initial conditions 
for the Big Bang can be solved unambiguously in the SM and in the i/MSM, and to find 
these initial conditions. To this end we consider in detail how the energy stored in Higgs- 
inflaton gets transferred to the SM and z^MSM degrees of freedom. This allows to make 
a refined estimate of the reheat temperature and to fix the concentrations of the singlet 
fermions before the hot stage. We show that the abundances of new particles are too small 
to infiuence the low temperature baryogenesis in the z^MSM studied in [|, IHl,^ and low 
temperature dark matter production worked out in ||l^, 14, p!s|].'^ 

The Lagrangian of the SM or of the z^MSM, which can be considered as the effective 
theories, can contain all sorts of higher dimensional operators, suppressed by the Planck 
mass. Therefore, we consider the infiuence of these operators on infiation and on production 
of singlet fermions of the z^MSM. We find that these operators are definitely not essential 
for baryogenesis and for dark matter production, if mass of the lightest sterile neutrino 
is below 100 keV. In other words, the conclusion that the production of WDM sterile 
neutrinos with mass in the keV region must be due to their mixing with active neutrinos 
is a robust consequence of the z^MSM. Since the presence of higher-dimensional operators 



^Perhaps, the observed accelerated expansion of the Universe should not necessarily be included in this 
list as it may be irrelevant for the early stages of the Universe evolution we are interested in this work. 
^See for a suggestion to use singlet fermion oscillations for leptogenesis. 

^For an original proposal of sterile neutrino as a dark matter candidate see |2^, |2^, earlier compu- 
tations of sterile dark matter abundance can be found in ^ ^t\ . 
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looks to be a generic phenomenon, we argue that the DM sterile neutrinos have to be 
lighter than few MeV in order not to overdose the Universe. 

The paper is organized as follows. In section § we review the mechanism of inflation 
based on the Higgs boson of the Standard Model, and determine the relevant interactions 
of the Higgs-inflaton with the other fields of the SM. In section |^ we analyse different 
processes reheating the Universe after inflation and estimate the maximal temperature 
Tjnax- In section ^ we discuss the approximate conservation laws in the SM and the i/MSM 
and define the operators Qa- Then we estimate the values of chemical potentials [xa 
generated by renormalizable interactions existing in the z^MSM and the SM. In section ^ 
we add to the theory higher dimensional operators and analyse their influence on Higgs- 
driven inflation and on the generated values of the chemical potentials. In section ^ we 
study the effects of CP-violation at the reheating stage. Section |^ contains conclusions. 



2. Higgs-driven inflation 



The inflationary model with the Higgs boson as the inflaton 
coupling with gravity to the action of the SM (or i/MSM) 



adds the non-minimal 



Sj = SsM + / d'^x 
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M2 



(2.1) 



Here S^m is the SM action, M is some mass parameter, which is nearly equal to the Planck 
mass in our case, R is the scalar curvature, $ is the Higgs doublet, and ^ is a constant 
fixed by the requirement of correct scale of the CMB fluctuations. Index "J" stands for 
the "Jordan frame" action. Action ( p.l| ) contains all possible terms of dimension 4 without 
higher derivatives.^ In this section we review shortly the inflation analysis of |21, ^] and 
introduce some formulas important for the study of the reheating period. 

The only part of the action relevant for inflation is the scalar sector. In the unitary 
gauge with $(x) = it has the form 



Sj 



d'^xyj—g 



M2 + ^/i2 ^^h^^'h 
R H 7^ 



„2\2 



(2.2) 



and Higgs vacuum expectation value is v = 246 GeV. Another part we analyse later on, 
while we always stick to the unitary gauge for simplicity. 

The conformal transformation. The simplest way to work with this action is to get 
rid of the non-minimal coupling to gravity by making the conformal transformation from 
the Jordan frame to the Einstein frame (see, e.g. [E^, |30|): 



Ml 



(2.3) 



*One could also add other dimension 4 terms like , R^^^R^'^ , etc., but they lead to terms with higher 
derivatives in the equations of motion and, therefore, lead to additional degrees of freedom, which should 
be dealt with in some special way. We do not consider such extensions here. 
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where Mp = I/VSttGat = 2.44 x 10^^ GeV is the reduced Planck mass. This transformation 
leads to a non-minimal kinetic term for the Higgs field. So, it is also convenient to replace 
h with new canonically normalised scalar field x by making use of 



dh V 04 
Finally, the action in the Einstein frame is 



/<f.V-S^-^/;+^-C(x)^ (2,5) 



where R is calculated using the metric g^y and the potential is rescaled with the conformal 
factor 

We will a bit ambiguously write potential U and scale factor Vt as functions of either h or 
X, which should not lead to misreadings, as far as h and x can be expressed one through 
another in a unique way. Figure || illustrates the connection between the Higgs field in the 
Jordan frame, /i, and the Higgs field in the Einstein frame, x- For ^ ^ 1, the solution of 
eq. ( p.4p can be approximated in two major regions,^ separated by 

_ [2Mp 

{h for /i < Xcr , 

JlMp\ogVt^{h) for X^,<h. '■^■^^ 



Namely, 



Note that for analysis of reheating we will need only the field values smaller than Mp / y^, 
where the logarithm in (|2.7| ) can be expanded in the following way 



Using relations ( |2.71 ) we can explicitly write the potential as (here we assume v <C Mp/^) 

ix* for X < Xc, , 

Uix):^{ w. A ^-^y (2,9) 



4|2 

Again, in the region interesting for reheating, the potential can be approximated by the 
quadratic potential 

U{x)^Y^'' for X„<x«yfMp, (2.10) 
where the "infiaton mass" uj is 

fXMp 

"-VaX- (2.11) 

Figure ^ shows schematically the potential (p.Sj). 

^Exact analytic solution exists, but is not really enlightening. 
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Figure 1: Dependence of x (the Einstein 
frame Higgs field) on h (the Jordan frame 
Higgs field), logarithmic scale. 



U(X) r 




Figure 2: Effective potential in the Einstein 
frame. The insert magnification is not to 
scale. 



Inflationary phase. The potential ( |2.9| ) is exponentially flat for large field values, and 
provides the slow roll inflation. Analysis of the inflation in the Einstein frame® can be 
performed in the standard way using the slow-roll approximation. The slow roll parameters 
(in notations of |j|]) are easier to express analytically as functions of the field h using ( |2.4 ) 
and ( p.6D , instead of the field 

Mj, fdU/dxY _Ml, fU' lY ^2.12) 



2 \ U J 2 \U x' , 
= Ml ^— ^ = Ml — , (2.13) 



where ' denotes derivative with respect to h. Slow roll ends at e ~ 1, which corresponds 



to the value /lend- The perturbation modes of WMAP |34] scale /c/oq = 0.002/Mpc left 
horizon when the field value equals /iwmap- The latter is determined by the number of 
inflation e-foldings, 

'^cnd r 

To generate proper amplitude of the density perturbations the potential should satisfy at 
/iWMAP the normalization condition 

U/e = 24tt^A'^mI ~ (0.0276Mp)'^ . (2.15) 

For usual quartic potential inflation this condition fixes the coupling constant A, while in our 
case this allows to find the value for for any given value of A. The infiationary predictions 
(see, e.g., Q) for the CMB spectrum parameters are then given by the expressions for 
spectral index Ug and tensor-to-scalar perturbation ratio r, 

ns = l-6e + 2r] , r = 16e , (2.16) 

also calculated at /iwmap- 



^The same results can be obtained in the Jordan frame jsl], |3^ . 
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In the case of the standard Higgs potential V{h) = ^(/i^ — f we get for the slow roll 
parameters |28| (in the limit h"^ > Mp/£^ S> v"^, ^ 1, exact expressions can be found in 



Hi), 

Inflation ends in at /lend - (4/3)i/'^Mp/i/f ~ 1.07Mp/^/^ (and Xend - 0.94Mp). The 
number of e-foldings is ( ^.141) 



^WMAP ^end , i„„ ^ + ^^cnd/^P 



+ log 



(2.18) 



leading to /iwmap ^ 9.UMp/^/^. Thus, the WMAP normalization pTTs] ) requires (for 



N = 59) 

C 47000\/A , (2.19) 

where A is the Higgs boson self coupling constant, taken at inflationary scale. Note, that 
we retained here the logarithmic term in (|2.18D , which was left out in [^]. This, together 
with WMAP5 value for normalization (|2.15| ), changed the numerical value in the relation 



( 2.19 ). This does not significantly change the spectral index and tensor to scalar ratio. 

The spectral index is ~ 1 — 8(4A^ + 9)/(4A^ + 3)^, and the tensor-to-scalar pertur- 
bation ratio is r ~ 192/(4A^ + 3)^. 

The number of e-foldings is fixed from the post-inflation history of the Universe 
described in section ^. We show there that the inflationary stage is followed by the matter 
dominated epoch, corresponding to oscillations of Higgs-inflaton with frequency to, defined 



in ( 2.11| ). The radiation dominated era starts at effective temperature Tr, given by ( 3.13| ) 



Then, the number of e-foldings is (see pjqj ) 



k , lOi^GeV ,^ C/1/^(XWMAP) 1, U^/Hx^ 



end ) 



~ 60.4- log -^-i log ^ . (2.20) 
Here the present Hubble parameter is Hq = 0.7/(3000 Mpc), ^7(xwmap) — ~Atf--, p{Tmax) 



is the energy density at the beginning of the hot stage, is in the range ( |3.1C| ). Then, 
we get 

iV ~ 59 , ~ 0.97 , r ~ 0.0034 . (2.21) 

The predicted values are well within one-sigma border of allowed region of parameter space, 
see figure |3[ 

Effective couplings in the inflationary domain. The inflation and reheating of the 
Universe occur at energy scales much larger, than the electroweak scale. This calls for the 
study of radiative corrections to the inflationary potential. A qualitative discussion of the 
influence of loop effects on inflation can be found in [^]. A number of explicit computations 



(giving in some cases conflicting results) has been reported recently ||36|, |38|, [40|. 



The conclusion of 38, ^ is that Higgs-driven inflation is a viable phenomenon in a 
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certain interval of Higgs masses mmin < < "^max• The values of m^i^ and rrimax found 
in 1 39] are: 



mt- 171.2 a^- 0.1176 

[126.1 + ^-^^ X 4.1 - p^^^ X 1.5] GeV , (2.22) 

[193.9 + H^l^^Il:! X 0.6 - ^l^HiZ^ X 0.1] GeV . (2.23) 

^ 2.1 0.002 ^ ^ ^ 



with a theoretical uncertainty (5theor = i2 GeV. The similar result for rrimin was found in 



an earlier paper |3S] and also in [40]. 

What is essential for the present 
study of reheating, is the magnitute of 
the coupling constants in the relevant 
energy domain ~ Mp/^. So, an appro- 
priate renormalization group running of 
the coupling constants should be taken 
into account. Specifically, one should 
use the Mp/^ scale value for the elec- 
troweak coupling constant 



a 



W 



43 



(2.24) 




Figure 3: The allowed WMAP+BAO+SN region for 
inflationary parameters (r, Ug), adopted from ]34| . 
The green box is our predictions supposing 59 e- 
foldings of inflation. Black and white dots are pre- 
dictions of usual chaotic inflation with A(/)^ and rn^cjp 
potentials, HZ is the Harrison- Zeldovich spectrum. 



and the corresponding numbers for the 
strong and U(l) gauge couplings. 

The dependence of the value of the 
scalar self-coupling A at the scale Mp / ^ 
on the Higgs boson mass is illustrated 
in figure § (see ref. [39[ for detailed de- 
scription). It can be seen, that out of 
the window (|]2|,|2]2|) for Higgs masses the inflationary scale Higgs self-interaction starts 
to behave badly at the energy scale of inflation. For large Higgs masses it becomes large 
and thus leads to strong coupling. For small Higgs masses it gets negative and leads to 
instability of the electroweak vacuum. The analysis of the present paper is not applicable 
very close to the boundaries of the allowed region. If the Higgs mass is approaching ( 2.22| ) 
or ( p. 231) one has to redo the analysis including higher order radiative corrections to the 
Higgs potential, what is beyond the scope of this paper. 

Yet another remark concerns applicability of perturbation theory at high momenta. 
As was found in [^l| and in [^, ^ , the perturbation theory, which is used throughout the 
calculations of reheating, is inapplicable for momenta above Mp/^ in a potential of the 
form ( |2.6| ). However, at reheating the energy is mostly contained in particles with momenta 
~ \Mp/^ (see below), which interact weakly. Thus, the details of the formulation of the 
theory for high momenta are irrelevant for present calculations. 



3. Reheating in Higgs-driven inflation 
3.1 Qualitative picture 
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The scalar potential for the Higgs 
field"^ in the Einstein frame exhibits three 
qualitatively different behaviours, leading 
to three stages of the Universe expansion. 
The first, inflationary stage, corresponding 
to the flat potential at x > Mp, has been 
already discussed in section ^. The second 
specific region of the scalar field values is 
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Figure 4: A at the scale Mp/^ depend- 
ing on the Higgs boson mass mn for rrit = 
169.1, 171.2, 173.3 GeV (from upper to lower 
graph). Varying e- foldings number and an error 
in the WMAP normalization measurement one in- 
troduces changes invisible on the graph. 



where the scalar potential is essentially 
quadratic, see (|2.10| ). The slow roll infla- 
tion terminates at x ~ Mp with the onset 
of the oscillations of the scalar field. Since 
the effective inflaton mass to is non-zero for 
these field values, the exponential expan- 
sion of the Universe is changed to the power 

low, corresponding to matter domination. The amplitude of the Higgs field during this stage 
is decreased due to expansion of the Universe and due to particle creation. At last, for 
X < Xcr we are certainly in the radiation-dominated epoch: the potential for the Higgs field 
(|2.9| ) does not contain any essential mass parameters and thus is scale-invariant; the scalar 
self-coupling and couplings of the Higgs field to the fields of the SM are relatively large, 
and lead to a rapid energy transfer from the coherent oscillations to relativistic particles. 
Assuming the instant conversion of the energy of coherent oscillations to relativistic degrees 
of freedom of the SM, we get a lower bound on reheating temperature Treh ^ 1-5 x 10^^ GeV 



(see m 



However, as we show in this section, creation of particles happens to be important 
even for x > X^r, and the reheat temperature is higher. In what follows we will be 
mostly interested in the very moment, when matter dominated expansion is replaced by 
the radiation dominated one. That is when energy in coherent oscillations of the scalar 
field is equal to energy collected by SM particles. We will characterise this moment by an 
effective temperature (the "r" stands for "radiation dominance"). It is determined by 
equating of the would-be thermally equilibrium energy of SM plasma described by to 
its actual energy. The real thermal equilibrium is achieved at somewhat lower temperature 

^reh ^ ^r- 

To determine T^, let us neglect first the effects of particle creation and consider the 
evolution of the Universe with the Higgs field in the interval ( p.ip . The Friedman equation 
reads: 



1 



3M|, 



a;2 ^2 



1 



(3.2) 



''We will use the names "Higgs boson" and "inflaton" interchangeably, depending on the context. 
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and leads to a matter dominated expansion regime 

a oc t^/^ , (3.3) 
X{t) = X{t) cos [uj{t - to)\ , (3.4) 

Here t is the physical time, a is the scale factor, H is the Hubble parameter, X{t) is the 
amplitude of the background inflaton field oscillations, to here gives the arbitrary phase 



of the oscillations, and uj is defined in ( 2.11 ). This solution is approximate. It is only 



reliable for H <^ uj, when the change of the scale factor is small during one oscillation. The 
amplitude reaches the critical value Xcr at the critical time 

t^t,, = -^ . 3.6 

UJ 



To determine the particle production, we will consider the solution ( |3.5D as an external 
background. This approximation breaks down when the energy of created relativistic 
particles is comparable with the energy of the scalar field (inflaton zero mode) 

Anf = = jXlx' . (3.7) 

This moment will give us the temperature we are interested in. 

To start with, we describe on the qualitative level various processes which occur during 
the reheating stage at t < tcr and single out the most important ones. Further, we analyse 
these processes in detail. 

The main mechanism draining energy from the inflaton zero mode is creation of the 
particles directly from coupling to the Higgs-inflaton. In the background approximation 



the inflaton field (3^) can be considered as an external source of all other fields. This 
source has the form of the varying-with-time masses of all the particles (this includes the 
propagating modes of the Higgs field itself). Only particles with large couplings to the 
Higgs field can be created effectively by this mechanism. These are the gauge bosons and 
top quark. Their masses in the region (2.8) are 



2 / X ff^ Mp\x\ 
m^x) = > (3-8) 

Mx)=m^^^si,nx. (3.9) 

Here /^tt = aw is the weak coupling constant, and yt = V2mt/v is top quark Yukawa. 
However, exactly due to large couplings they are heavy and are still non-relativistic. The 
production of such particles does not change the equation of state from the non-relativistic 
matter to radiation. That change happens eventually only due to creation of the relativistic 
secondary particles (such as light leptons or quarks) via decays or scatterings of the heavy 
particles. A competing (but slightly slower) process is the direct creation of relativistic 
Higgs excitations, which happens because the potential ( ^ ) is not exactly quadratic near 
the origin.* 



Sr 



The non-linearity of the potential (2.9) at large field values x ~ is relevant for particle creation 



only during a short period at the very early time, that may be neglected. 
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As a result, the generic picture of the reheating process is the following. As far as 
the inflaton "mass" uj is smaller than the gauge boson (^^) or top quark mass ( |3.9| ) for 
X ^ Xcr, creation of the gauge bosons or top quarks is possible only at the moments, when 
the inflaton field crosses zero (when xit) ^ A^cr)- During each zero crossing some gauge 
bosons and top quarks are created. At first, when the concentration of the created particles 



is small (occupation numbers nk <C 1), the creation rate is constant (see Appendix A. 2 ). 
At this stage the created W bosons are non-relativistic and decay into light SM fermions 
(which are relativistic). The decay rate, however, changes with time with the decreasing 
amplitude of the inflaton oscillations. The decay process sustains some quasi constant 



density of the created bosons (Appendix BJ). It stops when the decay rate becomes smaller 



than the production rate, which happens at the inflaton oscillation amplitude (|3.21 ). Up 



to this moment no significant energy transfer from the inflaton to radiation happens. Then 
the generation process accelerates, being enhanced by the stochastic parametric resonance 
(occupation numbers > 1, and the concentration of W bosons rises. The energy transfer 
into the light SM fermions proceeds now mainly via WW — > // annihilation, (see Appendix 
|B.2| ). This process rapidly transfers all the energy to radiation, resulting in transition from 
the matter domination expansion a oc t"^^^ to the radiation domination a oc at field 



amplitude only slightly smaller than ( 3.21 ). This should be considered as a refined upper 
bound^ on critical X. The lower bound is given by the slower energy transfer mechanism — 
the generation of the Higgs bosons on close-to-vicinity nonlinearities of the potential ( |2.g| ) 
(see paragraph ^ , which yields transition at the oscillation amplitude ( p.27| ) . We do not 
analyse the production of the top quarks in the present work, because their contribution 
is smaller than that of W bosons. This is because the parametric resonance enhancement 
is absent for fermions due to Pauli exclusion principle. 

To summarise, the matter-radiation transition happens when the amplitude of the 
inflaton oscillations is somewhere in the region 

1 /2 

The temperature Tr is estimated as follows, 

where g^: ~ 100 is the effective number of degrees of freedom of the SM. This gives for 



/ \ \l/4 

1.4 X lO^^Mp <Tr< 4.5 X 10"^ Mp , (3.12) 

\0.25 / 

or 

3.4 X 10^3 GeV < < l.lxlO^^GeV. (3.13) 

The coupling constant A here is taken at the inflationary scale. Its dependence on the 
physical Higgs mass is presented in flgure ^. Note, that at T = the particle distributions 
are not yet fully thermal. 



^More accurate analytical estimate of outcome of the scattering processes is hardly possible, since particle 
masses vary quite rapidly. 
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3.2 W boson production 

Let us start with description of boson production by the external oscihating source (^.41). 
We write the equation of motion with the mass given by ( |3.8D . Actually, with formula 
( |3.8| ) we ignore the exact behaviour of the mass in time intervals when x < ^cr- then the 
potential is quartic, so the zero mode evolution also deviates from (3.4). Hence, the jump in 
the derivative in m^r is "smoothed" . This introduces an additional cutoff in the spectrum 
of produced W bosons, which can be safely neglected. To simplify the computation, we also 
replace the vector boson by a scalar particle. We will take into account three polarizations 
of the vector boson in the final formula only. 

Evolution of the mode (/^k with conformal momentum k is governed by the equation 



(3.14) 



This equation can be solved in the adiabatic approximation except for the moments when 
myy{t) is close to zero. The reason is that mwit) is larger than the frequency of the 
background field lo except for small background x ^ -'^cr- The solution of equation ( 3.14| ) 



in the adiabatic approximation is (see Appendix A.l) 



l3/2 



•■ Jo hidt 



+ 



Pi 



hi Jq kgdt 



(3.15) 



where parameters a^, remain constant between the moments tj, corresponding to zero 
background field x{tj) = 0. 

In the vicinity of the moments x{tj) = t^is mass m\y{t) becomes small compared to 
the background (source) frequency u and particle creation can take place. Exact solution 

) allows for mai 

in the occupation numbers = |/3kP we have 



in this region (see Appendix A.l) allows for matching (^^,(3^^ to Ok'/'^k- ^'^^ change 



5nl 



l^kl 



+ 



2|i?k| 



l-|i?k 



2 + Y 



2|fik| 



nlinl + 1) 



(3.16) 



where R\r is a decreasing with k function defined in ( A. 15] ) (see also figure ||) , and 
the k-dependent phase ( |A.17 ). 



tot 



IS 



If the occupation numbers on the right hand side of eq. ( 3.16 ) are small, <C 1, 



it describes a simple particle creation with constant rate. This is certainly the case if the 
particles, generated at each zero crossing, decay before the next zero crossing, or scatter and 
change momentum k to some value where the generation coefficients in (3.16) are small. If 
not, eq. ( 3.16| ) describes a resonance like production of bosons. There exist several regimes, 
depending on the model. The usual ( "narrow" ) parametric resonance emerges when rij ^ 1 
and the phase OI^^/tt changes by an integer number between tj and tj+i. In this case one 
has a pronounced exponential behaviour for rik, located in narrow regions of momentum, 
which are defined by A^:^^^. = ^tot^ ~ ^tot — Another resonance situation is realized 
when 

(3.17) 



This regime is called "stochastic resonance" ||4j]. In this case the jumps at moments tj 
may be in either direction, but, on average, they also lead to exponential growth, though 



AC » ^ • 
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a slower one. The regime (|3.17D is realised for X > {\/g'^)Xcr (see (|A.17D , (|A.8D ), which is 



always true in our case. In this regime the particle generation proceeds for any momenta 
with not too small |-Rk| (otherwise the last term easily spoils growth of the occupation 
number), which implies low momenta. 

So, while the number of created particles is small, <C 1, the creation of the 
bosons from ( 3.16| ) proceeds in the linear regime and can be approximated as (see Appendix 



,^3.a3^^^2^2x, (3.18) 
dt 27r2 ' ^ ^ 

with numeric coefficient A ~ 0.0455. The created particles are essentially non-relativistic. 

For concentrations of other gauge bosons we have the obvious relations nw+ = n^-, 

nz = ny^+ 1 cos^ Q-w^ where By/ is the weak mixing angle. 

When the occupation number becomes larger, the production is enhanced due to the 

Bose statistics, and is becoming approximately exponential (Appendix |A 

d{(j?n\^ 



dt 



a-^2LoBnw , (3.19) 



where the numerical coefficient B ~ 0.045. These particles are also created with nonrela- 
tivistic momenta. Note again, that here the backreaction of the particles on the condensate 
is neglected. 

3.3 Transfer into relativistic particles 

Now let us analyse how decay and scattering of the W bosons influence their generation, 
described in the section 

At small concentrations of the bosons the main process is their decay. The (average) 
decay width of the SM gauge boson is 

r^0.8aw{mw), (3.20) 

The mass here is approximated as averaged over inflaton background field oscillations. 
Comparing ( |3.20| ) with the production rate ( 3.19| ), we find that for 



^ > nL K^ X,, ^ 40 (-^ ] X,, . (3.21) 



0.64 TT " \0.25 

the decay is a more rapid process and prevents the exponential regime to start. 



The energy transfer rate is then balanced by the linear production ( 3.18 ) and the decay 
rate (|3^) . In Appendix pT| we show that the energy transfer to the relativistic modes is 
negligible for this period. 

When the decay process becomes inefficient, the concentration grows, leading to en- 
hancement of the production approaching exponential behaviour ( 3.19| ). At the same time 



the main process responsible for the energy transfer to light particles becomes the annihila- 
tion of the W bosons, which is proportional to density squared. The relevant process is the 
annihilation into two fermions WW ff via t-channel fermion exchange. The estimate 
of the cross section for the process is 



" ^ [7=2) — — (<T ^ '° W ' ^'-''^ 
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where Ni = 3 is the number of lepton generations, A'^g = 2 + 1/4 is the effective number of 
quark generations (virtual t-quark contribution is suppressed by (miy/m-t)^ ~ 1/4), and 
A'^c = 3 is the number of colours. 

Then, the equality between generation and annihilation of the W bosons is reached at 

- ^■^'^ 

'^scatter — • [o.Zo) 

a 

The energy drain into the relativistic modes is 

d 



dt 



a 



V) = 2o^ • J(m2^) • anLtter • (3-24) 



The integral is saturated at late times and gives nearly immediate transfer of all the energy 
into the relativistic modes after the regime (|]2l|) finishes. 

One may note, that this approximation may overestimate W boson production, and it 
may actually proceed at a smaller rate, because of several reasons. First, the occupation 
number nk is not too large (|B.9|) , and resonance regime is not fully reached. Second, the 
exponent in ( 3.1S| ) is actually the upper limit. Careful analysis may reveal that the process 



is slower, so the transfer to the relativistic degrees of freedom happens later (at lower X) 



than (3.21), the latter should be considered as the upper bound on Xr- The lower bound 



is then given by a slower process of generation of the Higgs field excitations. 
3.4 Higgs production 

Another particles produced during the inflaton oscillations are inflaton excitations (Higgs 
particles). Really, mass of the excitations 6x in the background x is given by = f7"(x), 
which is, approximately 

f for x{t) > , 



where x(^) is given^'^ by (^!^). Thus, we need to solve the same equation (|3.14| ), but now 
with the mass ( 3.25| ). The details of the solution are given in Appendix j^. The produced 



particles are relativistic, with energy E ~ IVSXX, and with energy balance equation 
density 

d(a^p) o (j^ , , 

-V^ ~ a^—^ . (3.26) 

dt 27r3 ^ ^ 

This competes with the inflaton (zero mode) energy density pinf at 

This provides the lower bound on the moment of transition to the radiation dominated 
epoch Xr- 



^"Of course, this is a rough approximation for X 3> ^cr, because we should in principle solve the equation 
of motion for x in the exact potential. But at large X S> Xcr the time, system spent in the region x < ^cr , 



is small, and the background solution can be simply approximated by {i-4) 
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3.5 Fermion production from Higgs decay 



For completeness, let us also calculate the number of light fermions generated at the re- 
heating stage by the inflaton-Higgs field. The result here does not apply to heavy fermions 
(top quark). The latter abundance is of small interest: before thermalization it is in any 
way not larger than that of the gauge bosons, and after thermalization top quarks are 
generated rather fast. 

We analyse here the production of the light fermions by the Higgs condensate decay 
due to the Yukawa interactions. The latter are of the form ( |3.9| ) with some small Yukawa 
y instead of yt- The exact treatment would require the solution of Dirac equation with 
the time-dependent mass (p. 9]). To make an order-of- magnitude estimate, we will replace 
the source oc sin(u;t) by the simpler one, oc sin(a;t). Though the spectra of the produced 
fermions are different (the spectrum is monochromatic for the sinusoidal source), their 
total numbers are similar. After this substitution the fermion time-dependent mass term 
becomes: 



Mp 



X sin(a;t)'i/'V' • 



(3.28) 



In the lowest order of perturbation theory the rate of fermion production here is equivalent 
to the decay rate of the system of scalar particles with mass u, concentration n = luX'^/2 

and effective Yukawa constant y^ 



d 3 .coX^ y2 ^Mp\ 

— (a n^] = a — ^to — — — , 



(3.29) 



(and the same formula for the antiparticles n^). An elementary computation leads to the 
constant physical particle density during the matter dominated expansion. 



2iO^Mp 1 



167r ^/SX ^ 32^/27r 



Xi PS 80 



A 



0.25 



y^Xl . 



(3.30) 



It is convenient to compare with the entropy after the end of the matter domination 
stage 



47r2 4vr2 / 30A . ^ 



2.9 



0.25 / 



100 



0.25 ) 



3/4 



xt 



(3.31) 



where we adopted (|3.12|) for the range of Tf. The resulting abundance is in the range 



Sr 



A 



0.8 ^ 28 1 I 

0.25 7 V0.25y 



A \ 



5/4' 



(3.32) 



These results are used below for an estimate of primordial abundance of the sterile neutrinos 
in the i^MSM. 
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4. Initial conditions for the hot Big Bang 



In section ^ we found that when the amphtude of the Higgs-inflaton drops below 3.8Xcr 
the matter-dominated expansion of the Universe is changed to the radiation dominated 
behaviour, which can be characterised at this moment by the effective temperature — 
3 X 10^^ GeV. It is this moment which can be considered as a starting point for the standard 
hot Big Bang: the later evolution of the system can be followed with the use of the SM or 
the z^MSM Lagrangian and standard finite temperature equilibrium and non-equilibrium 
methods. As we discussed in the Introduction, to specify the system completely, one has 
also to determine the values of the chemical potentials corresponding to either exactly or 
approximately conserved quantum numbers. In this section we identify the most important 
operators and fix the chemical potentials for them. 

Let us start with the SM. It has got three anomaly-free exactly conserved quantum 
numbers Qa = — {L^ is the lepton number of the generation a and B is the bary- 
onic number). In addition to them, there are quite a number of approximately conserved 
different fermionic numbers, such as baryon number (broken by electroweak anomaly), 
asymmetry in the number of right-handed electrons and light quarks such as u and d (bro- 
ken by small Yukawa couplings), etc. In the standard inflationary logic one concludes that 
all the quantum numbers — eigenvalues of corresponding operators — are exponentially 
small at the end of inflation (at h ~ Mp) and thus can be put to zero. What concerns 
the charges Qa, they cannot be created in the process of reheating the Universe, anal- 
ysed above, simply because they are exactly conserved. As for the other charges, such as 
asymmetry in the light quark flavours, their generation can only occur due to CP-violating 



effects. Therefore it is suppressed by the Jarlskog determinant |45|, since the only source 
of CP-violation in the SM is related of the Kobayashi-Maskawa phase. Applying the argu- 
ment of refs. [46, ^ to this case one concludes that asymmetries in all CP-odd operators 
at the beginning of the Big Bang are at most on the level of W~^^. To summarize, in the 
SM all chemical potentials are negligibly small at the beginning of the Big Bang. At the 
same time, the CP-even operators (such as the abundance of fermions plus antifermions 
of a given type) equilibrate with the rate not smaller than ce^T, which exceeds the rate 
of the Universe expansion right after the beginning of the Big Bang. So, deviations from 
thermal equilibrium in the CP-even sector in the SM can be neglected as well. 

Let us turn now to i^MSM. The most general renormalizable Lagrangian containing 
the SM fields and three right-handed singlet fermions has the form: 

£ = £sM + Njid^j'^Nj - LaNj^ - ^ NfNj + h.c. , (4.1) 

where £sM is the Lagrangian of the SM, Uai are new Yukawa couplings, $ is the SM Higgs 
doublet, and = ^ij^'j- In the contrary to the SM, there are no exactly conserved global 
quantum numbers in the model. To analyse the approximately conserved currents, at least 



the orders of magnitude of various parameters, entering eq. (4.1), have to be fixed. In 
the z^MSM the Majorana masses of the singlet fermions are below the electroweak scale 
and, correspondingly, the new Yukawa coupling constants are smaller than those in the 
quark or charged lepton sector. Their values are constrained by cosmology, astrophysics 
and experiment. 
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The hghtest out of the 3 neutral leptons, A^i, plays the role of the dark matter particle. 
Its Yukawa couplings are bounded from above as 



^|F„i|2< 10-24 (4.2) 



from cosmological considerations |Q, Q related to DM production and from X-ray con- 
straints on the radiative width of the DM sterile neutrino []To| . 

The heavier nearly degenerate singlet fermions A''2,3 (their common mass is denoted by 
M) fix the pattern of neutrino masses and mixings and produce baryon asymmetry of the 
Universe. Their Yukawa couplings satisfy: 

where matm ~ 0.05 eV is the atmospheric neutrino mass difference, k = 1(2) for normal 
(inverted) neutrino mass hierarchy, = [y'^y]ii, and e = ^ < 1. If e ~ 1, i.e. for the 
case when the couplings of singlet fermions to active leptons are similar, F2 is at most 
~ 2 X 10"^^, corresponding to M ~ M^y. For the smallest possible value of parameter 
e ~ 7 X 10^^ (a lower limit is coming from the requirement of successful baryogenesis, see 
figure 10 of [^]) one gets an absolute upper bound on F2, 

F| < 3 X 10-^^ , (4.4) 

roughly coinciding with the electron Yukawa coupling. 

In the limit Fj ^ the sterile fermions completely decouple from the fields of the SM, 
and the i/MSM contains an infinite number of exactly conserved operators, corresponding to 
a number of singlet fermions with any given momentum. Right after inflation these numbers 
are exponentially small and can be put to zero. Then, these singlet fermions are created as 
described in section ^.^^ With the use of (|4^ , [4^ , |3^) we get that the abundance of DM 



sterile neutrinos A^i produced at the reheating stage is at most Ai = nAr-^/s~7x 10"^^, 
and the abundance of A''2,3 is at most 2 x 10"^. These numbers are too small to play any 
role in the subsequent evolution of the Universe. These were the constraints on CP-even 
operators, the CP-asymmetries in left-right helicities are suppressed much stronger as the 
CP violating amplitudes must contain at least two extra powers of Yukawas. 

To summarize, the initial condition for the Big Bang in the z^MSM can be described 
by the density matrix 

/5(0) =/)sM® |0)(0| , (4.5) 

where psM = exp(— /3?t!sm)i P = ^/T, is the equilibrium SM density matrix at a tem- 
perature T with zero chemical potentials, and |0) is the vacuum state for sterile neutrinos. 
The physical meaning of eq. ( |4.5| ) is clear — it describes a system with no sterile neutrinos, 
while all the SM particles are in thermal equilibrium. It is this expression which was used 
for computation of DM abundance and for computation of baryon asymmetry in the z^MSM 
in refs. H, |l|, |1,|,|U. 



"We work in the basis in which A/12 = M13 = 0, A/23 = A/32 = A/, A/i = A/n < A/, A/22 ~ M33 < M. 

^^Yet another mechanism for production of singlet fermions is the decays of Z and W bosons to sterile 
neutrino and left-handed lepton. Since the rate of this reaction is suppressed not only by the square of the 
same Yukawa coupling but also by an extra gauge constant, we expect it to be subdominant. 
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5. Higher dimensional operators 



In the first part of the paper we analysed the inflation and reheating in the model with 
action ( p.ll ). However, one may expect that there are corrections to this action, suppressed 
by some large energy scale (the Planck mass) . We consider the following higher dimensional 
operators as an addition to the SM or z/MSM Lagrangian: 

<5£nr = ^^^^^N'^N - ^($^^)' + ^Ll^^^L, + • • • + h.c. . (5.1) 

The natural value of all dimensionless coupling constants is about one. 
The following questions arise: 

1. Do these operators spoil the picture of inflation discussed above? 

2. Does the reheating change? 

3. Can the singlet fermions be created due to these operators in substantial amounts? 

In this section we analyse these issues. Along the lines of consideration in section |2|, 
we rewrite first these operators in the Einstein frame. For Higgs part 5Cmr this yields a 
modification of the potential in (p.2|), ( |2.6| ). The transformation rule for the Higgs- fermion 
interaction is readily obtained, if we also make the conformal transformation of all the 
fermionic fields ip 

iIj^'^ = 17-3/2^ . (5.2) 

The kinetic part for the fermions is conformally invariant, while the Yukawa part of the 
action 

'S'j,Yukawa = / <fx^/^ Y{h)'l(jlp , (5.3) 



changes. Here Y(h) describes the generalised Yukawa interaction providing the fermion ip 
with mass. It is not important for the present discussion whether this is the Majorana or 
the Dirac fermion. The corresponding Einstein frame term is 

Si5,Yukawa = J cl^Xy^ ^(^^(^^ ^ " ^^'^^ 

Specifically, for the Dirac mass this yields the Einstein-frame terms like 



which we already used in deriving ( |3.8| ), (p.S|). Note, this mass rescaling is similar to that 
of massive gauge bosons. 

For the Majorana higher dimensional term in (|5.1|) we get 



<5/:i.,NR,Majorana = ^^^'^^ " (5-6) 

Note, that at the reheating stage, (Mp/^ <C x ^ Mp/^/^), we have ~ 1, so the only 
change is the field substitution in accordance with (|2.8D. 
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5.1 Contributions of the higher dimensional operators to inflation 



Clearly, the main effect of the higher dimensional operators is expected when the Jordan 
field h is large, of order Mp. Indeed, at this scale the higher order Higgs field operators 
may spoil the fiatness of the Einstein frame potential (2.9), and the fermion Majorana 



mass terms can also give sizeable radiative corrections to the potential. This can change 
the infiationary properties of the model. In particular, inflation may turn to be impossible 
if the sign of the slope of the potential is changed in the inflationary region. Or, the pre- 
dictions of the CMB spectral index and tensor-to-scalar ratio may leave the experimentally 
admitted region. However, the constraints on higher dimensional operators, imposed by 
the requirements that the inflation is not spoiled, turn out to be rather weak. The reason 
is that the inflationary potential is only essential at sufficiently small values of the Higgs 
field, h < hwMAP ~ lOMp/v^. These values are well below the Planck mass, so that 
non-ronormalizable contributions are well suppressed. We analyse below in some detail the 
contribution of operators (|5.1|) to infiationary potential. 



5.1.1 Higgs operators 

Let us analyse the following higher order terms added to the Higgs potential: 



n = 6, 8, . . . . As far as the operators are suppressed by the Plank mass, their effect is mostly 
important at high values of h. At this scale we have-*^^ dx/dh ~ \/^Mp/h, ~ ^"^h^/Mp 
The contributions to the slow roll parameters at > 60 are 



2(n - an ( h X""^ 

The main contribution comes from the lowest order power term, . Thus, for the change 
of the parameters at the normalized-to-WMAP value of the field we have 

5e ~ 1.7 X IQ-^al , (5.10) 
6r] ~ 0.00506 , (5.11) 

for A = 0.25. This implies the change in the spectral index 5ns — O.Olag — O.OOOlog and in 
the tensor to scalar ratio 5r ~ O.OOOSog. 

To keep the spectral index within la bounds 0.94 < < 0.98 (at small r), see 
figure ^, the coefficient for the dimension six operator in the Higgs potential should be 
lag I ^ 3. Hence, no significant contributions are expected from the higher order operator 
with natural values of the coefficients of order one. 



13 



At the end of inflation the values of dx/dh and Q, are different, but it only slightly changes the WMAP 



value for ^, as far as the contribution from (5.7) are more suppressed for lower h 
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5.1.2 Yukawa terms 



Let us now analyse the effect of the Yukawa terms for the sterile neutrinos on the inflation. 
They come from the fermionic loop contributions to the effective potential for the Higgs 



field. According to (5.4) and (5.6) the mass term for the right handed neutrinos in the 



Einstein frame has the form 



Mi ^ 



2^{h) 2MpQ.{h) 



where M/ is the usual Majorana mass term for the sterile neutrinos in z^MSM. For large h 
the first term is suppressed, but the second one (dimension 5 operator) provides the mass, 
growing with the field. The latter could change significantly the effective potential for the 
Higgs field. Indeed, the term yHLN induces, in the Einstein frame, the usual Dirac lepton 
mass, the growth of which stops at /i ~ Mp/y^, being suppressed by Q{h). This dimension 
5 contribution to the mass yields the following contribution to the Higgs effective potential 
in the inflation region 



m 



4 (h) m%{h) P^h^ 



= _^,„,^.__^,„,_. (5.12) 

For high enough /3 this changes the sign of the derivative of U{h) at some h, which would 
make the inflation impossible or limit its duration. Let us calculate the slope of the 
potential 

U Ux' " h^^ V 8^^^^^ / X' ' 

where ' means derivative with respect to h and we neglected the derivative of the logarithm. 
It is required that U'{h) > for at least 60 e- foldings of inflation, i.e. for all h satisfying 
h\/^/Mp < 10. The logarithm here, accounting at least for the inflationary epoch, 1 < 
h\/^/Mp < 10, is about log(lOO) ~ 5. This implies the constraint 

/3Mog(/iV^^V) ^-,.-6 



87r2Af 



< 10" 



leading to 



/?^<47(^). (5.14) 



For smaller Higgs masses the bound is stronger, but always much larger than one. We 
see, that for rather large value of the dimensionless constant in front of dimension-5 mass 
operator for the right-handed neutrinos, the inflation is not spoilt. It is also straightforward 
to check by exact calculation of the spectral index that constraints from the WMAP on 
lead to essentially the same bound on fi. 



^''Of course, one could imagine starting inflation exactly from the top of the potential, where the derivative 
of U{h) is zero. However, this corresponds to a highly tuned situation, keeping in mind that the change of 
the derivative is due to interplay between tree level term and radiative corrections to the effective potential. 
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5.2 Sterile neutrino production 



As we have seen in section 3.5, during preheating the sterile neutrinos are produced very 



slowly by the renormalizable dimension 4 operators. Let us estimate the contribution of 



the dimension-5 operators to sterile neutrino production during and after preheating. 
We will separately analyse the production in the thermal bath after reheating by the 
annihilation process hh — > A^A^, and during reheating by the decay of the inflaton-Higgs 
condensate. 

5.2.1 Thermal production 

Let us start with the study of neutrino production in the primordial plasma. In this 
section we consider the neutrino production after reheating of the Universe, at T ~ Tj. 
(the higher temperatures are more essential due to the suppression of the relevant operators 
by the Planck mass). Here the electroweak symmetry is restored, and the production of 
sterile neutrinos goes through annihilation of the Higgs bosons (4 degrees of freedom, 
corresponding to the unbroken phase of the SM) due to coupling (|5.1|). 
The cross section of this process is (neglecting the neutrino mass) 



In the absence of other sources of neutrino production, the interaction (|5.6| ) contributes to 
the r.h.s. of the Boltzmann equation for sterile neutrino density 

^(o^nAf) = a^Aahh^NNnl . (5.16) 

Here we took into account the annihilation of all four modes, re/j stands for the density of 
each scalar degree of freedom. 

This equation can be easily integrated accounting for the fact that at T < the 
Universe is at the radiation dominated stage: 

a^Vi, nh = , (5.17) 



H=- = \ —^ = — , s = g^ , (5.18) 

a V 90 Mp 2t ' ^ 90 ' ^ ^ 



The solution of equation (5.16) gives the neutrino density-to-entropy ratio. 



135C'(3)^/l0 (3^ 

where g^ = 106.75 is the total number of degrees of freedom of the SM. We also set the 
initial abundance to zero. We will compute it in the next chapter, dealing with sterile 
neutrino production during preheating. 



^^We suppose that at the temperature Tr the Universe is aheady in a fully thermalised state. Though 
this is not exactly the case, we expect that this assumption can only overestimate the number of produced 
neutrinos, since the non-equilibrium Higgs spectra are more enhanced in the infrared region in comparison 
with the thermal one. 
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Putting the numbers, we get at low temperatures T <^Tr: 



Ajv = 1.5 X IQ-^P^ 



Mp 



(5.19) 



With the use of this relation we can answer the question whether the primordial thermal 
production of the lightest practically stable sterile neutrinos can substantially contribute 
to the DM abundance. The neutrino-to-entropy ratio remains intact, so at the current 
moment we have 



So 



where = 0.87 x 10^^*^ is the baryon-to-entropy ratio and nB,o is the present baryon 
number density. Therefore, we can write for the sterile neutrino abundance 0,^ 



N 



0,B f^N A 



N 



where mN and rup are the sterile neutrino and proton masses, respectively, 0,b 
and r^DA/ = 0.23 are the baryon and DM abundances |34]. Hence, 



N 



Mn P^Tr 



Qdm 27keV Mp 



0.046 



(5.20) 



So, for the maximal allowed (3 ( ^.14 ) and for the reheat temperature in the range 



( 3.13 ), we conclude that the neutrino mass, required to provide the proper DM abundance, 
should be in the range 



Mn 



0.25 



13 



0.25 



1/4 



MeV ^ 42 MeV 



(5.21) 



If the dimension 5 operator is present with the "natural" coefficient /3 ~ 1, then the mass 
of the (long living) sterile neutrino should not exceed 



'0 25\ -^/'^ 
Mn < 600 MeV ( J 



(5.22) 



in order not to overproduce the DM (the upper bound from (3.13) is used for the estimate). 



5.2.2 Production during preheating 

In this section we consider the neutrino production in the early Universe right after inflation 
got terminated. The production during this period happens due to the effective interaction 
with the Einstein frame field x 



C 



JL 



IxlA^^'iV + h.c. 



(5.23) 



We will see that this mechanism produces more sterile neutrinos, than the thermal pro- 
duction discussed above. 

To find particle production due to the time-dependent fermion mass one has to study 
the Dirac equation following from ( ^.23 ). However, to simplify the discussion, we proceed 
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as in section 3^ and replace \x\ by x- Then the rate of fermion production coincides, to 
the lowest order in Yukawa coupling, with the decay rate of a collection of scalar particles 
with certain mass and number density. An analysis performed in Appendix ^ shows that 
the number density of produced fermions is not affected by this replacement, though the 
spectrum changes. 

We can write the Boltzmann equation as 

where we replaced the oscillating source by an ensemble of free scalar particles with the 
number density = ^X'^ and particle decay width F^^attv = 

Then with background (p.4|), ( p.5| ) one gets (the early-time contribution is negligible): 

Dividing this by the entropy ( 3.31| ) we get 



^ /32 ^30\^/'l [xTr Mp 



32Tr^ \XgJ 32^2 j; 

i.s.:o-(2?)^. (...) 



This is larger than contribution from thermal generation ( 5.1S| ), for the reheating temper- 
ature in the range ( |3.13| ). 

Proceeding analogously to the previous subsection we get 



Qdm 2.2xl08keVV X J T, 
So, for the maximal allowed /? given by ( p. 14 ) and for the reheat temperature in the range 



( 3.13 ), we conclude that the neutrino mass, required to provide proper DM abundance. 



should be in the range 



Mn = 65 keV ^ 210 ( ^ j keV . (5.26) 



If /5 ~ 1, that is what is naturally expected, we have 

Miv < 3 (^^) MeV . 

Let us summarize the results obtained above. If a more complete, than the i^MSM, 
theory leads to higher-order non-renormalizable operators characterised by a "natural" 
constant /3 ~ 1, then the mass of the DM sterile neutrino must not exceed few MeV.^^ 
Otherwise, it will be produced in amounts enough to overdose the Universe. Sterile neu- 
trinos, produced at reheating, can only play the role of CDM, since their mass must exceed 



'For smaller /3 this limit scales as Mjv oc l//3^ 
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65keV. This requirement comes from the inflationary upper hmit on /? ( 5.14| ). Finally, 
if the sterile neutrino has a mass in 0(10) keV region and thus plays a role of WDM 
candidate, the thermal primordial production, discussed in this section, plays no role. 

The higher dimensional operators, of course, produce also heavier singlet fermions of 
the z^MSM. In section ^ we analyse whether this has any influence on the low-temperature 
baryogenesis due to singlet fermion oscillations. 



6. Higher dimensional operators and baryon asymmetry 

It is shown in section |5| that the abundance of DM sterile neutrino, created at reheating 
due to higher dimensional interactions, cannot exceed 

An ~ 10"^ . (6.1) 

Hence, the "primordial" (related to inflation) creation of DM is not effective for light 
sterile neutrinos, and may play a role only if > 65keV. Interestingly, this number is 
only somewhat larger than an upper limit on the mass of DM sterile neutrino produced 
resonantly due to lepton asymmetry generated in the z^MSM |l^, |l9|. Mi < 50keV. In 
other words, the initial condition (^]^) is certainly valid for sufficiently light singlet fermion 
(mass below 65 keV), which could play a role of WDM. 

Other singlet fermions can be produced due to the same type of higher dimensional 
interactions (note that for the abundance computation the magnitude of the Majorana 



mass plays no role), and their abundance is bounded from above by (6.1). This number 
is much smaller than one, meaning that the heavier singlet fermions are practically absent 
at the beginning of the hot Big Bang. Still, in order to proof that the density matrix ( [4.5| ) 
can be used as an initial condition, one must show that the CP asymmetries in distribution 
of singlet fermions do not exceed the baryon asymmetry of the Universe ~ 10"^''. 

It is not difficult to see that this is indeed the case. To this end consider the most gen- 
eral form of the leptonic part of the Lagrangian, taking into account the higher dimensional 
operators of dimensionality 5 as well: 

Ccp = j^^^^N'jNj + i^L^M^Lp - y^j L^Nj^ + 5a/3La^/3^ + h.c. , (6.2) 

where Ej^ are the right-handed charged leptons. To get an amplitude of CP-violating 
effects, one may consider the imaginary parts of re-parametrisation invariant products of 
Yukawa couplings, which can be considered as a generalization of Jarlskog invariant for the 



Kobayashi-Maskawa quark mixing to this case (see also |46, ^). These invariants can be 
written as traces in flavour space of the products of fa/3, Vai and gap (no contraction 
between Greek and Latin indexes). 

The fermion production due to Higgs oscillations and Higgs scattering appears first to 
the second order in these couplings. As we have seen in section ^ the leading effect comes 
from the first term in ( |6.2| ). Clearly, there is no CP-violation in this order. To the fourth 
order in coupling constants the CP-violating effects appear through the CP-violating trace 
Tr[y"'"/3y/]. Since from the flatness of potential |/| ^ /? ^ 6, and because yai are strongly 
bounded from above by ( [4.2| , [4.4| ), there is a suppression of the asymmetry at least by 10 



orders of magnitude in comparison with (|6.1|) . Going to higher orders makes the situation 
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even worse. To conclude, the initial conditions for the Big Bang are correctly described 
by eq. (4.5), even if higher dimensional operators are included in the i^MSM, and thus 
the baryogenesis is a low-temperature phenomenon, having nothing to do with inflation or 
Planck scale physics. 



7. Conclusions 

In this paper we analysed in detail the evolution of the Universe in the scenario where 
the Higgs boson of the SM plays a role of the inflaton. The history of the Universe can 
be divided into three stages. The first one is inflation. Here the non-minimal coupling of 
Higgs to gravity makes the effective scalar potential flat, the Universe expands exponentially 
and the necessary spectrum of perturbation is generated. This stage finishes roughly at 
h ~ Mp. During the second stage the Universe expands as under matter domination. 
The Higgs field oscillates in the nearly quadratic potential for Mp/^ < h < Mp, and the 
particle production is not effective. When h reaches the critical value h ~ Mp/^ the energy 
stored in Higgs zero mode is transferred rapidly in other degrees of the SM, producing the 
hot Big Bang with temperature ~ 10^^ GeV. After this time the Universe is dominated 
by radiation. 

We have shown that at the onset of the radiation dominated epoch the densities of all 
CP-odd operators in the SM can be put to zero and demonstrated the for the case of the 
i/MSM the concentrations of the singlet fermions are negligible at T^. 

We also considered an extension of the SM and I'MSM adding to them higher dimen- 
sional operators suppressed by the Planck scale. We analysed the constraints on these 
operators coming from the condition to have successful inflation. We demonstrated that 
the concentrations of the singlet fermions at can be safely put to zero, provided the 
mass of DM sterile neutrino does not exceed 100 keV. This means that in this case the pro- 
duction of baryon asymmetry and of dark matter must occur at small temperatures (about 
and below the electroweak scale) by essentially the same mechanism, as was described in 
1 18, The properties of singlet fermions can be almost unambiguously fixed by different 



cosmological considerations ||T^, . 

We found that the presence of higher-dimensional operators provides a new mecha- 
nism for primordial production of DM sterile neutrino. This mechanism is effective in 
models with sufficiently heavy sterile neutrinos, M^r ^ 100 keV. No presently available 
astrophysical constraints (in particular, those associated with X-rays) can exclude this 
possibility, since production occurs even if DM sterile neutrino Yukawa couplings are iden- 
tically equal to zero. However, if the solution of the short scale difficulties of the CDM 
scenario 49, 50, 52, ^] is to be given by the WDM, this region of the parameter 



space should be discarded. At the same time, there are no reasons to expect that these 
operators are suppressed by the scale exceeding the Planck one (i.e. it is unlikely that 
/3 < 1). Therefore, the models with DM sterile neutrinos heavier than 3MeV are generally 
disfavoured due to problems with dark matter overproduction. These arguments provide 
an extra justification of sub MeV mass of the lightest singlet fermion within the z^MSM. 



^'^See, however, ref. where it is shown that heavy sterile neutrinos could be WDM for other types of 
production mechanisms. 
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Note added 



Some time after our paper was posted at arXiv the article devoted to the same subject 



appeared. Most of the conclusions of |56| are similar to ours. In particular, the authors 
of used a similar formalism to analyze the transfer of energy from the Higgs field 
oscillations to gauge bosons. A detailed analysis of differences and similarities of these 
works goes beyond the scope of the present paper. We would like just to mention that 
some (not inessential for physical consequences) differences in numerics are presumably 
due to the fact that annihilation of created gauge bosons was not accounted for in p6| , 
leading to the different rate of transfer of energy to relativistic particles at later stages of 
reheating. 
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A. W boson production 
A.l Semiclassic approach 

The discussion in this section closely follows [^]. At first approximation a creation of 
W bosons can be regarded as a creation of particles with mass ( |3.8D varying with the 
amplitude of the Higgs field ( |3.4| ) . This approximation breaks when important amount of 
energy is transferred from the inflaton zero mode ( |3.7| ). If at this moment the energy is 
in the relativistic modes this corresponds to the moment of transition to the radiation 
dominated expansion. 

To solve the equation ( |3.14| ) we rescale the variables by 

= a'/'</'k . (A.l) 

This leads to the equation for an oscillator with varying frequency 

^k + A:g(t)^k = 0, (A.2) 

kl{t) = — + ml,{t) + A, (A.3) 

where A = — | (^) — |^ is always small and can be neglected. The initial conditions, 
corresponding to vacuum oscillations, are 



^k(i) = . (A.4) 



Equation ( A.2 ) can be solved in the adiabatic approximation when /cq ^o- At k = this 



condition is equivalent (up to a change of the scale factor a, which is negligible in our case) 
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to rhyy <C mi^. This is true for 



\t - tj\ > 



1/3 



1 



2g'^MpXuJ 



41/3K ' 



(A.5) 



where tj are moments when inflaton crosses zero, so that Tn?{tj) = 0, and K is the natural 
scahng parameter defined in ( [A.12 ). In these regions adiabatic solution is 



+ 



(A.6) 



where parameters a^, remain constant within tj-\ < t < tj. At the moments tj the 
coefficients get changed by the Bogolubov transformation 



al+^e-< 




(A.7) 



where -Rk and -Dk are the "reflection" and "transition" coefficients for each interval < 
t < tj (they obey the equality + l-D^P = 1), and 



9^ 



kodt . 



(Ai 



To find the coefficients i?k and Z^k for each interval we need to solve exactly eq. ( [A .21) 
in the vicinities of the moments when rriw ~ (and where the adiabatic approximation 
is inapplicable) and match this solution with ( [A.6| ) in the intermediate region of the field 
amplitude. Obviously, exact solution is impossible, but we can approximate m^(t) near 
zero as const - {t — tj\. This approximation to the potential is good enough at 



, V6 
\t-tj\ < — 

CO 



The regions ( |A.5D and ( |A.9D intersect for 

X > 



A 



48\/67ravy 



Xr 



(A.9) 



(A.IO) 



which covers all the possibly interesting reheating period. At the same time, one can 
estimate, that the effects of "smoothing" of the \t — tj\ are also insignificant up to approxi- 
mately Xcv Thus, we can match the solution of the linear equation and adiabatic solution. 
To solve the "linearised" equation we rescale the variables as 



with 



K 



6^' 



1/3 



g2 MpX{tj 



2^/6 e 



1/3 



(A.ll) 



(A.12) 
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Then, for small r eq. ( A. 2 ) takes the form 



dr2 



+ {k^ + IrD^k = 



(A.13) 



which can be readily solved analytically in terms of the Airy functions. Matching the 
solution with the asymptotic form ( |A.6| ) at t = tj one gets (see Appendix ^ for details) 



2^(i^+.) i [Ai^(-^^) Bi(-K^) - M{-k') Bi-(-K^)] 

(Bi(-K2) + i Ai(-K2))(Bi'(-K2) + i Ai'(-K2)) 
- Ai'(-K2) Ai(-K2) _ Bi(-K2) Bi'(-K2) 



4 y 



(Bi(-K2) + i Ai(-K2))(Bi'(-K2) + i Ai'(-K2)) 



Thus, we can calculate the occupation number = at the moment tj 



n 



\ , ^ ,o ni + 2\ 1 + ni\ ni 



a] 
^tot 



i^kl 



cos(6l^ot) 



-20^ - 2 ( + ^ ) + arga^k - arg/3^ 



(A.14) 
(A.15) 



(A.16) 
(A.17) 
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Figure 5: Particle creation coefBcients for (3.16), and the effective resonance exponent parameter 
(A. 23). The right plot is in the log-log scale, to show and k,~^ behaviour at large momenta. 

The total particle number density for tj < t < tj+i is given by 



{27ra)" 



n{t, < t < t,+i) = / 7T^< . (A.18) 



A. 2 Non-resonance production 

In this case, we estimate the production of the particles simply by the first term in ( |A.l(j ) 
(only one degree of freedom of W boson is accounted for) 



dt IT J (27r)3 |L>kp 27r3 
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with the numeric coefficient A = 0.0455. The spectrum of the created particles is presented 
in figure ^. It has cutoff at |k|/a ~ K and power law tail |k|^^. Thus, at X > Xcr the 

produced particles are nonrelativistic, as far as K < (rra^). So, to study the transition 
to the radiation domination one should analyse further the W boson decays into relativistic 
particles, which is done in the next subsection. 

Let us also note, that taking into account three polarizations of the vector bosons and 
their different types is made as follows 

3 

ny/+ = n^- =3n , nz = TT^n , (A. 20) 

cos"^ Ow 

where dw is the weak mixing angle, and n^yi, nz are densities for and Z-bosons, 
respectively. 



A. 3 Stochastic resonance 



If the occupation numbers rik exceed one, then the first term in ( |A.16 ) can be neglected 



and the last two terms, proportional to nk, yield the resonance — an exponentially rapid 
particle creation. If ?ik ^ 1 and the shift in the phase, AO, is large, then we can neglect 



the first term in (A. 16) and write approximately 

n^+^ - n^e^'^'^kavg . (A.21) 

The average growth exponent /^kavg is obtained by averaging of the exponent over the 
random phase 

The integral can be found exactly, leading to 

^..M = ^MMi , (A.23) 

which is presented in figure |5|. 

If we neglect the expansion of the Universe, then the time derivative of the total particle 
number can be estimated as 



/d'^k / |k| \ 

^^nk2a.^..,^— j~2..i?n, (A.24) 



dn f d'^k 
'di 



where the numerical coefficient B ~ /Uavg(O) ~ 0.045. The exact expression depends on 
the exact spectrum of the generated particles, and can be omitted at our level of precision. 
Note, that the same equation, describing the exponential creation, is true without any 
change for n\Y± and nz- The difference appears in the pre-exponential behaviour only. 

Typical momentum of the produced particles is again K, so they are non-relativistic, 
and analysis of their conversion into light relativistic particle is needed. It proceeds via 



annihilation (scattering) and is studied in detail in Appendix B.2. 
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B. W boson decays and scatterings: energy transfer into relativistic par- 
ticles 



B.l W boson decays 

Now let us analyse whether some processes, like decay or scattering of the W bosons, may 



destroy the resonance picture described in Appendix |A.1| . These processes may destroy 
the resonance behaviour by taking the bosons out of the resonance region in two ways: 
either by transferring the energy to other particles, or by changing the boson momenta 
and taking it out of the resonance region. 

Changing momentum could be expected in WW — > WW scatterings. However, as the 
typical W boson momentum is smaller than their mass K < (mw), one can not achieve 
in scatterings momenta larger than K, and this process (though rather effective) can be 
safely neglected. 

So, the two remaining processes are decay and annihilation of the W bosons, which 
transfer the energy to the relativistic (light) particles, and, depending on their rate, may 
also prevent the development of parametric resonance. 

We start with the analysis of the decay process of the gauge bosons created at the 
moment tj. The (average) decay width of the SM gauge boson is given by (|]20|). We also 
estimate the W boson mass as the averaged value over inflaton oscillations, 

/ 2 \ 9^ Mp{\x\) ^ Y y .-R T ^ 

= ^ ^ = ^Oiw^crX , (B.l) 

where {\x\) = — ■ If the decay is faster than the exponent of the stochastic parametric 



resonance 2loB (see ( 3.19 )), then the parametric resonance never settles, and creation is 



dominated by the first term in ( [A. 161) , see section [A. 2 . The inequality F > 2ujB leads to 



( |3.21[ ). So, for the period before ( |3.21| ) the production happens only due to the first term 
in ( p.l6| ). Let us check, that the energy in the decay products of the W bosons remains 
small for this period. 



During the time period ( 3.21 ), when the decay is fast, the creation of the particles is 



non-resonant ( p.l8|) . We can write the approximate Boltzmann equation for this period 

The solution to this equation in the semi-stationary regime, corresponding to vanishing 
time derivative in the left-hand side, is 

ndecay ^ • (B.3) 

The semi-stationary approximation \h\, 3Hn <^ Tn, is valid for X < ^^^(0.8aw)'^(,'^Xcr ~ 
0.7 X W^Xcr, that is always after the end of inflation. 

We can also check, that the occupation numbers are really much smaller than one 
and we are in the non-resonant regime ( p. 18 ). As far as the typical physical momenta of 



the W bosons are of the order of K we have 

ndecaj _ \/Act a a/a n na f ^ 
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which is much smaller than one for all interesting X. 

The energy during this stage is converted to the SM particles produced in W boson 
decays. They are light and relativistic, and typical energy transferred to them in each decay 



is of the order y (m^), as far as the W bosons are non-relativistic. Thus, we can write the 
Boltzmann equation in the expanding Universe for the energy density in the relativistic 
SM particles as 

3 \ r—r-AujK^ 



where the coefficient ( 6 H — I accounts both for the different number and masses of 

\^ cos-^ 8w J 

created W^- and Z-bosons (see eq. ( |A.20D ). The solution is saturated by late time and 
reads 

3 \ I—^AojK^ 6 



This reaches the inflaton energy density p ~ pinf (3.7) at 



- 6 + o . 2 4,1 ~ 5.8Xcr , (B.7) 



that is much later than the end of the non-resonant creation period ( p. 21 ). We conclude, 



that the energy drain by W boson decays is irrelevant during the non-resonant inflaton 
decay. 

B.2 W bosons annihilation 

Self scattering of the W bosons, like WW — > WW is of little interest for us, as far as it does 
not take the bosons out of the stochastic resonance zone (the bosons are non-relativistic, 
so after scattering they retain their small momenta). 

It is easy to check, that for the W boson number density (|B.3| ) saturated by the boson 
decays (discussed in Appendix |B.lj ), during the period ( 3.211) of non-resonant inflaton 



decay the annihilation to fermions is negligible, cm\^^^^ < Fndecay (relation formally holds 
for X > O.lXcr). 

The scattering proceeds much more actively at larger particle densities, so the rel- 
evant production mechanism is given by the stochastic resonance ( p. 19 ). Thus, we can 



approximate the effective Boltzmann equation for the W boson particle number as 

^ {a^nw) = {2Bujnw - crn^) , (B.8) 

where nw = n'^+ = riyy- . Its approximate solution ( |3.23| ) is obtained, again, by setting the 
derivative in the left-hand side to zero. This is true for d{a?' Uscatter) / dt <C a^2i?u;nscatter) 
that is for X < 45^Xcr « 4x lO^Xcr, while we are interested in much smaller X. We should 
check of course, that the particle density is not to small, to allow for stochastic resonance 
to work. Again, for the typical occupation number (up to some numerical factor) we get 

nk ~ = - 3.4 . (B.9) 
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This is larger, than one, thus we may hope that exponential creation is already a reasonable 
approximation for ( [A.16 ). The energy drain is obtain if we recall that the W bosons are 



nonrelativistic, so each scattering provides the energy transfer of 2-\/(m^), leading to 



(3.24). The solution of ( 3.24| ) is also saturated at late times, so (if the initial time is small) 
we have 



73 



A 



65 vr ^Jaw \0.25 
One can see, that it reaches inflaton energy density p ~ pinf (| 

2 2 „ 4 
1663 ^3 ^3 



H) at 



X 



2 2-1 

653 vr 3 A3 



-Xr: 



IWXc 



(B.IO) 



(B.ll) 



This is earlier, than the end of the non-resonant production region ( 3.21 ). Taking into 
account Z bosons makes this process even more active. This means, that after the moment 
(3.21), the parametric resonance starts, and due to higher concentration of the gauge 
bosons, the energy is rapidly transferred into relativistic SM fermions via gauge boson 
annihilation. So, we expect that the transfer of the energy to the relativistic modes via the 



gauge bosons completes by approximately (|3.2lD . 



C. Non-resonant Higgs production on the nonlinearities of the potential 
for smah x 



One needs to analyse the production of particles by ( 3.14]) , but with the mass 



(t) 



to 



for X cos [ujt] > , 



3XX^ cos^ [ujt] for X cos [cot] < 



(C.l) 



One way to find the generation by this source is to use the method described in section 
A.l. The adiabatic approximation holds while the mass does not change, and close to 



the moments tj the problem can be solved after replacing the cosine with the quadratic 
function. Alternatively, if the number of generated particles is small, a simpler perturbative 
approach can be used. We use the perturbative approach here. 

Let us first neglect the expansion of the Universe during several oscillations of the 
inflaton field. In this case the number of particle of the mass uj generated by the quadratic 
potential 



(t) 



UJ 



-•int 



is given by 



nk (t) 



■ UJ 



2ikot 



where /cq 



k"^ +uj'^. The integral is equal to (for the moment of time around t 



ti = ^l, 

/ = 1, 2, . . . ; the integral value changes while the inflaton field crosses zero, but the exact 
form is not important) 

2 / 21^^^ 



nk (ti) 



1 



sm 



16 sin2 
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where 



L 



3AX2 w 



1 



■ sm 



47r ( ^ + 1 ) e 



1 



2 kQ [kQ + to 
and the parameter e is defined from the equation 

sin [27re 



ko — LO 



sm 



4^ ( ^ - 1 1 6 

UJ 



UJ 



At large times, t 3> a; ^, using the equahty hm, 



3XX 

sin^ x-t 



(C.2) 



(C.3) 



t— >oo 



TTXrt 



5 (x) one gets 



Thus, integration over momenta gives a convergent sum 
n (t) 9X^X^ ^ 1 



257r3 

1=2 



1 



1 



1 



sin [47r (/ + 1) e] 



/ - 1 



sin [47r - 1) e" 



(C.4) 



The sum is saturated for Anne ~ 1. This imphes the typical energy of the produced 
particles, 

UJ 1 ; 

E V3AX , 

47re 2 

which is larger, than uj, so the particles are relativistic. 

Using formulas from Appendix |^ we get the following estimates for the production 

rate 

n (t) Attuj'^ uj 



n 



t 157r3 



and for the energy flux 



P 



p{€^^ UJ 



5 

27r"3 ■ 



Reintroducing the expansion of the Universe in the usual way by changing fi 
p — > we have the number and energy densities at late time 

47ru;^ UJ t 



n 



157r3 4 V3AX(t) ' 
3 uj^ 



11 27r3 



t . 



d{a^n) 
a^dt ■ 



(C.5) 
(C.6) 



D. Tunnelling through a — |x| barrier 

The solution of eq. ( [A.13 ) is given by the Airy functions for negative and positive times: 



^^{t < 0) = A_ Ai(r - K^) + B_ Bi(T - k^) , 
^k(T >0) = A+ Ai(-T - K^) + B+ Bi(-T - k^) . 



(D.l) 
(D.2) 
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The coefficients should be determined by the matching conditions at r = 
<^k(0-) = 'Pk(0+) , <(0-) = <(0+) . 



(D.3) 



It is comfortable, however, firstly to match the coefficients A±, B± with a = a^e '^k, 
P = /3^e+*^k, a' = a{'^^e~*^k, and (3' = /3j^'*'"^e*^k from the expansion ( |A.6| ). The asymptotic 
expansions of the Airy functions are 



Ai(-a;) 
Bi(-x) 



^ , sin I -x"'" + , 



i—rr cos (-X"'^ + , 
7rxi/4 \ 3 4 



-.3/2 , ^ 



(D.4) 
(D.5) 



Then, the solution matched with ( |A.6 ) at large r (for matching one should use in ( [A.6 ) 
only linear part of the mass, mw(t) ~ const - {t — tj\) 



^k(r < 0) 



^>k(r > 0) 



«e-(|«^+f)+^e^(§«^+f) 



/vr 



vr 



+ i 



a'e'^s"- ' 4^ + /3'e 



Bi(r - 

Ai(r - K^) , 
Bi(-r- + 
Ai(-r - K^) 



(D.6) 



(D.7) 



Using the linear relations between a, f3, a', (5' from condition at zero ( p. 3D , and the 
definition of i?ki from ( A.7), which is 



, 1 -R?l 



(Di 
(D.c 



one gets the expressions ( |A.14 ), ( A.15 ). 



E. Useful sums 

In Appendix ^ we obtain the following sums 



5i = E ^ l/l^ • sin [4vr (/ + 1) e] - ^ sin [4^ (/ - 1) e])' 

°° 1 / r / I 1 \^ 

1=2 ^ ^ 



(E.l) 
(E.2) 



where e is a small dimensionless parameter, e <C 1. These sums are saturated at I ~ l/47re, 
hence to get the leading order results in e one can replace these sums with corresponding 
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integrals by introducing a new variable u ss u = Anle, du = Aire. Then to the leading order 
in e one arrives at 

Si = r (2ucosu-2smuf = — (47re)^ , (E.3) 
Jo 15 

^2 = r e'^^ {2ucosu-2smuf = (Ane)^ . (E.4) 
Jo u 



F. Nonresonant peirticle production with h and \h\ sources 

Here we compare production of fermions by h and \h\ sources and conclude, that the 
corresponding production rates are the same, though the spectra differ. 

As far as the number of the created particles is small, nk <C 1, we can use the pertur- 
bation theory. Then, the perturbation 



leads to the number density at the moment t 

ft 



rik {t)= [ m {t') e^'^'Ut' [ m {t") e'''^'" di!' . 
Jo Jo 

and total particle number 

f d^k 

Here we calculate the particle density side by side for two different sources 



m{t) = msm{u!t) , (F.l) 
m{t) = \msm{Lot)\ , (F.2) 



and massless fields so /cq = fc- One can check that for n = 0, 1, 



2'r(n+2) 

/ dt sm (cot) e"'*"* = -—75 ^ • e* cos — • e* , (F.3) 

27r(w+l) 



Hence for the full {n + l)th period 

27r(n+l) 



I dt sm (ujt) e^*''* = —75 ^ • e' cos — • sm — e* , (F.5) 

7212 Ak^ -u^ V ^ / V / 



27r 



27r(rt+l) 

dt |sin {ut)\ e"'^' = ■ cos^ ( — ) • , (F.6) 





. ( kn 


(") 


• sm — 




V ^ 








1 • - , 
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and summing over + 1 periods one arrives at 



dt sm{ujt) e^'^'^o* = — ^ ^ -e - -cos — • sni — 



4A;2 - a;2 V 7 V ^ 7 sin (^) 

(F.7) 

I ■ / 2ifcnt . 2,rfc(iV+l) 2/«^^\ 

|sm(a;t)| e^*''"* = — ^ g • e • cos^ — 



4fc2 - u;2 V w y sin (^) 

(F. 

Hence the number of produced particles for + 1 periods of oscillations is 

• 2 / 2-Kk(N+l) 

27r(A^ + l)\ m2 /■ „ IQuJ^k^ o f kT^\ 2 ^'"^ I ^ 

n\t = ^ = — T / dk cos^ — • sin'' ' ' ^ 



a; / 2tt^ J {4k^ - uj^f \ io J \ lo J sin^ (2^) 

(F. 

k{N+l)\ 



. 2 / 2Trk(N 

n t = ^ = — / d/c ^ cos^ — ^„ ^ . (F.IO) 



Assuming that a tiny amount of particles is produced per each period one can turn to 
continuous variable T in these expressions. To obtain the particle production rate we are 
interested in linear in T contribution. It comes from poles in the integrands. Having this 
in mind and making use of the relation 

lim ^1^^ = 5(a) , (F.ll) 
one proceeds with calculations. For the source (F.l) one makes use of the identity 



while for the source ([F.2|): 



. 2 f 27rk{N + l) \ . 2 

sm = sm (kl ) . 

\ UJ J 

Then for the first source one obtains 

m2 f dk k"^ u? sin2 {Tk) 



Tfl P 



and with help of ( |F.11| ) 



, , m? f dk 7rA;2 ^ ( . uj 
2 



2^2 j (^^..^2 4 V 2 



This is the monochromatic spectrum. The number of produced particles is 

n{t)=t 



327r 
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For the second source ( p^.2| ) one has 



nit) 



r dkk'^ cos2 (^) sin^ (TA:) 



Here the double-pole at /c = 01/2 is cancelled by the double-zero from cosine squared. But 
there are a lot of poles due to sine in the denominator. In this case the useful variant of 

( IfH] ) is 



sin {tk) ijj^ 



It gives for the spectra 



lim „ . = -^b (k) . 

TTtsin^ 7r2 

dk k"^ — 



vn? iJ^ 



21 f) 

Finally, integrating over momenta and summing up the series 



n=oo 

E 

n=0 



n 



TT 

T 



one gets 



n{t) =t 



327r 



the same answer as for the first source ( [F.l| ). 
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